New High-Tc Half-Heusler Ferromagnets NiMnZ (Z = Si, P, Ge, As) 
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Based on the first principle calculation, we propose a new class of high-Tc half-heusler ferromagnets 
NiMnZ (Z — Si, P, Ge, As). The structural and magnetic properties are investigated through 
the calculation of the electronic structure, phase stability, equilibrium lattice constant, magnetic 
exchange interaction Jij and Curie temperature Tc. It is found that all alloys show half-metallicity 
and ferromagnetism at temperatures much higher than room temperature in a wide range of lattice 
expansion (compression). At the equilibrium lattice constant, Tc of 715K, 840K, 875K and 1050K 
are predicted by Monte Carlo simulation for NiMnP, NiMnAs, NiMnGe and NiMnSi, respectively. 
Following these results, these alloys are strongly expected to be promising candidates for spintronic 
applications. 

PACS numbers: 



I. INTRODUCTION 

The prediction and development of new magnetic ma- 
terials and half metalhc ferromagnets are necessary for 
exploiting the great potentials of spintronics in the near 
future. From the view point of industrial application, 
since the electronic structure behaves like metals w.r.t. 
the electrons of one spin direction and like semiconduc- 
tors w.r.t. the spins in the opposite direction (also called 
half-metallicity), leading to 100% spin polarization at 
Fermi level Ep, the half- metallic fcrromagnet with Curie 
temperature (Tc) higher than room temperature is con- 
sidered as a key material for spintronic applications. 

The family of half-Heusler alloys NiMnZ (Z is an 
sp— valence element) that have Clb structure as illus- 
trated in Fig. [1] are also greatly expected to be a fer- 
romagnetic half-metal because the interaction between 
the d orbitals of Ni and Mn (with suitable Z elements 
such as Sb) can form a gap between the bonding and 
anti-bonding states flj. The half-metallicity is first pre- 
dicted in the half-Hcusler alloy NiMnSb by de Groot 
and collaborators^] in 1983. Experimentally, besides the 
fact that the strong spin polarization has been verified in 
some compoundsfJU, the ferromagnetism in the half- 
Heusler NiMnZ at temperatures ranging between 500K 
and 730K has been known for Z== Co, Pd, Ft and Sb@]. 
Theoretically, the attempts to replace Ni (or Mn) of 
NiMnSb by other metals such as Fe, Co, etc. give a half- 
metallic picture in the band structure of some alloys [J 0] 
such as FeMnSb, CoMnSb and FeCrSb, etc. 

In this paper, based on the first principle calculation, 
we propose a new class of high-Tc ferromagnets NiMnZ 
(Z = Si, P, Ge, As) by placing the atoms lighter than 
Sn such as Si, P, Ge and As into Z sites of NiMnZ al- 
loy. For this aim, we have carried out the calculation as 
follows. First, we perform the total energy calculation 
to estimate the equilibrium lattice constants (ELC) and 



TABLE I: Structural configurations of NiMnZ. 
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configuration stability of Clf, structure by employing the 
ultrasoft pseudo-potential (UPP) method as the imple- 
ment in STATE-Senri code (thanks to Prof Y. Morikawa, 
Osaka Univ.). To evaluate ELC, we employ the general- 
ized gradient approximation (GGA) within UPP as well 
as both of the muffin-tin (MTA) and atomic sphere ap- 
proximations (ASA) within framework of LSDA (thanks 
to AKaiKKR package coded by Prof. H. Akai (Osaka 
Univ.)) [Sec. 2]. Using the obtained ELCs, we next cal- 
culate the density of states (DOS) [Sec. 3], magnetic 
exchange interaction between magnetic sites [Sec. 4]. 
Finally, using Jy as the input data we calculate Tc by 
employing three approaches: mean field approximation 
(MFA), random phase approximation (RPA) and Monte 
Carlo simulation (MC) [Sec. 5]. The effect of the lattice 
expansion/compression on the half-metallicity, magnetic 
exchange interaction and Curie temperature is also in- 
vestigated. The structural and magnetic properties of 
the prototypical half-heusler alloy NiMnSb are also cal- 
culated for comparison. 



II. PHASE STABILITY AND EQUILIBRIUM 
LATTICE CONSTANT 



We have performed the total energy calculation us- 
ing GGA within the framework of UPP method to es- 
timate the stability of three different configurations of 
NiMnZ: a, (3 and 7 phases. These phases are defined 
in Table |T1 The total energy differences between these 
phases show that the a phase is the most stable con- 
figuration and the (3 phase is the least stable one. For 
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FIG. 1: Clb structure: the unit cell contains Y at 4a(0, 0, 0), 
X located in the octahedral coordinated pocket 4c(|, j, j), Z 
at the body center 46(|, |, ^) and a vacant site at 4d(|, f , f )• 



TABLE 11: Equilibrium lattice constants calculated by em- 
ploying the MTA (aA/r) and ASA (a^s) within KKR-LSDA, 
and by GGA within UPP (aup)- Experimental value taken 
from Ref . I VII is given in the parenthesis. 



Alloy 


aMT (a.u.) 


OAS (au) 


aup (au) 


NiMnSi 


10.4696 


10.3291 


10.2612 


NiMnP 


10.5878 


10.4376 


10.2782 


NiMnGe 


10.7447 


10.5388 


10.4448 


NiMnAs 


10.8276 


10.6526 


10.5778 


NiMnSb 


11.5085 


11.3020 (11.2098) 11.0900 
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FIG. 2: DOS of NiMnSi. Left panel displays the partial DOS 
of 3p— Si (multiplied by 5), 3d— Ni and 3d— Mn, and total DOS 
at a = a AS- Right panel shows the total DOS within LSDA 
(KKR) and GGA (UPP) at three values of equilibrium lattice 
constant a — aup, a as and aMT- 



example, total energy differences for NiMnSi are given 
as: E{I3) - E{a) = 2.813 eV and Eij) - E{a) = 1.11 eV. 
Following this result, we next study ELC and magnetic 
properties of the alloys NiMnZ in the a phase. 

ELCs obtained by three approaches: GGA within 
UPP, MTA and ASA within KKR-LSDA are shown by 
Tab. HH It should be noted that GGA often gives ELC 
larger than LDA. However, the LDA within KKR ap- 
proach in our calculation produces ELC larger than GGA 
combined in UPP method. For the same compound, 
omt is the largest and ajjp is the smallest. Generally, 
the biggest deviation of umt (or aup) is approximately 
smaller than 3% in comparison with a as- In particu- 
lar, for the prototypical alloy NiMnSb, the experimental 
lattice constant [aexp = 11.2098 a.u.) is found to be in 
the range oi aup < a < a as and is close to the lattice 
constant obtained by ASA (a^s — 11.3020 a.u.) with 
a deviation smaller than 1%. Also, one can expect that 
the real lattice constant of these proposed alloys might 
be arranged in the range from amt to aup. 



III. DENSITY OF STATES 
A. Half-Metallicity 

DOSs at ELCs are shown in Fig. [21 In this paper, only 
DOS of NiMnSi will be presented as a typical case. The 



discussion of the half-metallicity in the remaining alloys 
will be presented elsewhere in more details. Moreover, 
one can also find a detailed explanation about the origin 
of the gap in various half-heusler alloys in Ref. IVII How- 
ever, it should be noted that although the DOS of NiMnZ 
with Z= As, Ge and Si were also considered in Ref. IVII 
the authors simply applied the experimental lattice con- 
stant of NiMnSb, that should be much larger than the 
real lattice constants (10% larger than aAs for NiMnSi, 
for example), for all of the alloys in question, resulting in 
the disappearance of the half-metallic character. On the 
contrary, using equilibrium lattice constants, we obtain 
the half-metallic character in the band structure of these 
alloys. The partial DOS (left panel) and total DOS (right 
panel) are calculated by LSDA and GGA at a as, clmt 
and ajjp- As seen from the left panel of Fig. [21 while 
most Ni-3d states concentrate at lower energies (in va- 
lence band) and are fully occupied, the 3d states of Mn 
are distributed to both lower and higher energies and 
are strongly polarized. The 3d— majority spin states of 
Mn are located at lower energies and almost occupied. 
They hybridize with Ni-3d states (and Si-3p states) to 
form a band at Ep . The minority spin states of Mn shift 
to higher energies and are unoccupied, leading to a gap 
(where Ep falls into) being formed in the minority band. 
Hence, the alloys are the half-metal. The calculated to- 
tal magnetic moment Mt of NiMnSi is exactly i^p in 
consistence with "rule of 18" Mt — Nt — 18, where Nt is 
the total number of valence electrons (iVj (NiMnSi) = 21). 
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Similarly, we have also obtained the half-metallic behav- 
ior and the total magnetic moment of 3fiB for NiMnGe, 
4fiB for NiMnP, NiMnAs and NiMnSb at a = aup, a as 
and qmt (except NiMnGe and NiMnSb at a = umt- See 
also Table Uni). 

B. Influence of Lattice Expansion on DOS 

To study the effect of the lattice expansion (compres- 
sion), we draw DOS of NiMnSi at different lattice con- 
stants in the right panel of Fig. [21 DOSs calculated at 
ELCs aijp,aAs and Qmt are plotted to illustrate the 
change of Ep according to the increase of lattice con- 
stant a. The figure in the top panel corresponds to DOS 
at the smallest a, and the figure at the bottom illustrates 
DOS at the biggest a. In addition, in order to compare 
with LSDA, DOS calculated by GGA (dashed hue) is 
also shown. The minority states of Mn-Sd electrons in 
GGA calculation shift toward the higher energies while 
the majority band is mostly kept unchanged, making the 
gap of the minority band wider than the LSDA. As the 
lattice constant a increases(decreases), Ep shifts toward 
the valence band (conduction band) and can fall into 
the valence (conduction) band, hence the half-metallicity 
might be destroyed if a exceeds the threshold value as in 
the cases of NiMnSb and NiMnGe at qmt- It is found 
that the half-metallicity of alloys NiMnSi, NiMnP and 
NiMnAs is preserved in the wider range of the expansion 
(compression) of the lattice cells. 

IV. MAGNETIC EXCHANGE INTERACTION 

To investigate the ferromagnetism in NiMnZ alloys, we 
calculate the magnetic exchange interaction Jij at three 
values of ELC. Once Jij is obtained, the statistical meth- 
ods such as MEA, RPA and MC are employed to calculate 
Tc- The exchange interaction Jij between two impurities 
at ith and jth sites, which are embedded in the ferro- 
magnetic medium, is efficiently calculated by utilizing 
the magnetic force theorem. The frozen potential ap- 
proximation [T^ is employed and Liechtenstein formula 
is used for evaluation of Jij 

J,, = — Im J dETT{{tl - t^)Gl^iE)(tl - tj)4(i?)}, 

(1) 

T / 1 

where t ■ are the site dependent scattering matrices for 

T / I 

spin direction up (t) or down (J,), and G-j are the KKR 
Green's hmctions connecting the ith and jth sites. The 
total energy change due to infinitesimal rotations of two 
magnetic moments at the ith and jth sites is calculated 
by using the magnetic force theorem, and the effective 
exchange interaction Jij is calculated via mapping the 
total energy change onto the classical Heisenberg model 
H = —Tii^jJijeiGj, where denotes a unit vector par- 
allel to the magnetic moment at the ith site. 
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FIG. 3: Mn-Mn effective exchange coupling constant Jij vs. 
distance d in units of lattice constant of NiMnZ (Z = Si, P, 
Ge, As, Sb). Jij is calculated at lattice constant a = aAS- 



A. Distance Dependence 

Fig- El displays Jij of NiMnZ at a = a as- In this figure, 
the influence of the sp-element on the magnetic exchange 
interaction is obvious. In alloys with iVj = 21 (NiMnSi 
and NiMnGe) , the exchange interaction is strong and fer- 
romagnetic for the 1st and the 2nd nearest neighbor pairs 
but anti-ferromagnetic for the 3rd and the Ath pairs. At 
farther distances, the interaction becomes weaker but fer- 
romagnetic. In alloys having Nt — 22 (NiMnP, NiMnAs 
and NiMnSb), the exchange interaction is ferromagnetic 
up to the Ath neighbor pairs but becomes weaker and 
anti-ferromagnetic at farther distances. In half-heusler 
alloys NiMnZ, J02 of alloys having a bigger Nt is stronger 
whereas Jqi is weaker. A similar dependence can also be 
seen in CMR manganitesjl^. 



B. Configuration Dependence 

The mapping on the Heisenberg model can produce dif- 
ferent sets of parameters in the calculation of Jy of differ- 
ent magnetic configurations (say, ferromagnetic and an- 
tiferromagnetic phases) for CMR manganites 13]. Thus, 
the evaluation of Tc based on sets can produce the 
different results for the same alloy. To consider the ac- 
curacy of this approximation in the calculation for our 
proposed half-heusler alloys, we now compare the elec- 
tronic structure and magnetic exchange interaction of 
two phases: ferromagnetic and disordered local moment 
(DLM) phases. The DLM phase is simulated by the con- 
figuration in which half of the local magnetic moments 
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FIG. 4: DOS of DLM phase in NiMnSi. 



point to the up direction and the others point to the 
down direction. DOS and Jij of DLM phase of NiMnSi 
are calculated in KKR-CPA-LSDA manner. As seen from 
Fig. [2] and Fig.gl DOS of DLM phase is quite similar to 
the ferromagnetic phase with the majority spin states of 
Mn concentrating at lower energies and minority spins 
of Mn mostly distributing in the higher energy region. 
Consequently, DLM phase exhibits a behavior of Jij sim- 
ilar to the ferromagnetic phase with a small change in 
strength (Fig. O. Jij of DLM phase is stronger than 
that of ferromagnetic phase for the odd (1st, Srd and 
5th) nearest neighbor pairs and weaker for the even (2nd, 
4th and 6th) nearest neighbor pairs. This differences 
is ignorable for the 1st and 2nd nearest neighbor pairs 
([J™ - J°^'^]/J™ is very small). Note that the mag- 
netic exchange coupling interaction in NiMnZ is short- 
ranged and very strong at very nearest distances and the 
dominant contributions into come from the 1st and 
2nd nearest neighbors; therefore, this change of Jij does 
not affect Tc so much. 



C. Effect of Lattice Expansion 

Now let us consider the effect of lattice expansion on 
the exchange interaction Jij . As shown in Fig. [21 the ex- 
pansion lattice causes the shift of Ep toward the valence 
band, hence it affects the exchange interaction between 
atoms on lattice sites. Fig. |6] illustrates the influence of 
the expansion lattice constant on the exchange interac- 
tion between Mn atoms of the fee lattice. NiMnSi and 
NiMnSb are chosen here as the typical alloys correspond- 



FIG. 5: Mn-Mn effective exchange coupling constant Jij vs. 
distance d in units of lattice constant of NiMnSi. Jij is cal- 
culated at lattice constant a = uas for the ferromagnetic 
(triangle) and disordered local moment (solid circle) states. 



ing to Nt — 21 and Nt = 22, respectively. of alloys 
with Nt = 21 (NiMnSi) and Nt = 22 (NiMnSb) are cal- 
culated at some values of the lattice constant a. For all 
cases in question, Joj {1 < j < 5) decreases with increas- 
ing a. The change of Jqj with higher j is very small and 
does not affect Tc so much. However, while Jqi of alloys 
with Nt — 21 decreases and its change is ignorable, Jqi 
of alloys with Nt — 22 remarkably increases with a. This 
leads to the different behavior of Tc when a varies: Tc of 
alloys having 21 valence electrons decreases whereas that 
of alloys with 22 valence electrons increases with a lattice 
expansion. Such a variation of Tc is shown explicitly in 
Table HTH 



V. CURIE TEMPERATURE 

After calculating Jy , the calculation of Tc is carried 
out by performing three statistical approaches: the mean 
field approximation, the random phase approximation 
and the Monte Carlo simulation. Tc is defined as 

ksTc^'^^^lY^Jo, (2) 

j 

in MFA, where fcs is Boltzmann constant, and by[l, [§| 

(fc^ri^^^)-l^ A^[J(0)_J(q)]-\ (3) 
q 
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FIG. 6: Jij VS. lattice expansion of half-heusler alloys with 
Nt = 21 (left figure) and Nt = 22 (right figure). The triangles, 
squares, stars and circles denote Jij calculated at aup, a as, 
aexp and umt, respectively. 



in RPA. Here, J(q) is the Fourier transform of the ex- 
change parameter Jij. 

To perform MC simulation, the Metropolis 
algorithm 1J| is applied to calculate the thermal 
average of the magnetization M and its powers. Then, 
the Binder's cumulant crossing methodflj] is employed 
and the fourth order cumulant C/4 defined by 



(4) 



is calculated as a function of temperature for different 
cell sizes (14 x 14 x 14, 16 x 16 x 16, and 18 x 18 x 18 
conventional fee cells) to find the universal fixed-point at 

In Tab. Illli we compare the Curie temperature Tc ob- 
tained by the Monte Carlo simulations with the results of 
MFA and RPA. It should be emphasized that the T*^"^ 
of NiMnSb calculated at ELCs as well as at the experi- 
mental value is in very good agreement with the exper- 
imental Tc = 730KH. It can also be seen that MFA 
overestimates Tc even at 100% magnetic atoms, whereas 
RPA underestimates Tc- Thus, as usual, the MC sim- 
ulation, which can be considered as the most accurate 
method, gives the Tc between values obtained by RPA 
and MFA at the same lattice constant. Also, the experi- 
mental Tc of NiMnSb lies between the values calculated 
by RPA and MC and is close to T^^^ . Hence, we will 
hereafter refer to the T*^*^ only. As expected, T*^'^ of all 
alloys is much higher than room temperature and ranges 
from 715K to 1050K (at qas)- The lowest Tc corresponds 
to NiMnP with Tc = 715K. NiMnSi has the highest Tc 
(1050K at a as) and the next highest ones are NiMnGe 
(875K), NiMnAs (840K) and NiMnSb (730K at ac^p). 
Due to the sp-element, the variation of Tc w.r.t. the ex- 
pansion of lattice cells is different according to Nt . Alloys 
with Nt = 21 have Tc decreasing with a (in the range 
where the half-metallicity is preserved), whereas alloys 



TABLE III: Tc calculated by MFA (Tj^^^^), RPA (T^^^-^) and 
MC {T^"^) at a = aMTjCiAS, aup and experimental lattice 
constant a^xp- The experimental value[^ of Tc is given in the 
parenthesis. 'HM' means half-metal. 



Alloy 




a 


rpMFA 


rpRPA 


rpMC 


Physical 








(K) 


(K) 


(K) nature 






aup 


1327 


1027 


1076 


HM 


NiMnSi 


21 


aAS 


1318 


1007 


1050 


HM 






aMT 


1198 


892 


968 


HM 






aup 


1178 


876 


942 


HM 


NiMnGe 


21 


a-AS 


1149 


834 


875 


*HM 






a.MT 


711 


561 




Metal 






aup 


746 


606 


625 


HM 


NiMnP 


22 


aAS 


863 


699 


715 


HM 






ClMT 


885 


723 


762 


HM 






aup 


944 


786 


816 


HM 


NiMnAs 


22 


aAS 


973 


807 


840 


HM 






a,MT 


962 


802 


867 HM 






aup 


877 


705 


736 


HM 


NiMnSb 


22 


Oiexp 


898 


716 


745 


HM 












(730) (HM) 






aAS 


928 


728 


748 


HM 






aMT 


839 


665 




Metal 



Although the LSDA calculation shows the 
half-metallicity, the results obtained by GGA present a 
metallic behavior in the bands of both spin channels. 



with Nt = 22 have Tc increasing. Accordingly, the half- 
metallicity in NiMnGe is considerably sensitive to the in- 
crease of the lattice constant due to Ep being very close 
to the valence band edge. In contrast to LDA calculation 
which gets the half-metallicity at a = a as, GGA cal- 
culation shows that the half-metallicity of NiMnGe can 
be destroyed if a > a as- However, similar to NiMnSb 
the real lattice constant of NiMnGe is expected to be 
smaller than qas and then the half-metallicity might be 
preserved. In short, the half-metallicity in the alloys pro- 
posed (except NiMnGe) can be preserved in a wide range 
of lattice expansion. 



VI. CONCLUSION 

In summary, we have investigated the structural and 
magnetic properties of the half-heusler alloys NiMnZ (Z 
— Si, P, Ge, As, Sb). The configuration stability is inves- 
tigated. The equilibrium lattice constants are predicted 
by three approximations within the ultrasoft pseudo- 
potential method and KKR-LSDA as well. Using the 
equilibrium lattice constants, the magnetic exchange in- 
teraction is calculated and then the calculation of Tc is 
performed by employing MFA, RPA and Monte Carlo 
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simulation. The role of sp-elements as well as the influ- 
ence of the lattice expansion (compression) on the half- 
metallicity and Tc are also investigated. With the results 
obtained, we propose new half-metallic half-heuslcr alloys 
NiMnZ (Z = Si, P, Ge, As) which have much higher 
than room temperature in a wide range of the lattice ex- 
pansion (compression). In particular, the present study 
shows that NiMnSi has very high Tc and hence is one of 
the most promising materials for spintronic application 
and we expect that it can be synthesized soon. In addi- 
tion, the high temperature half-metallicity and ferromag- 
netism can also be found in some Ni-based half-heusler 
alloys with Mn being replaced by Crp^. 
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